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Abstract. We report the discovery of fourteen new ZZ Cetis with the 4.1 m 
Southern Astrophysical Research telescope, at Cerro Pachon, in Chile. The can- 
didates were selected from the SDSS (Sloan Digital Sky Survey) DA white dwarf 
stars with T e ff obtained from the optical spectra fit, inside the ZZ Ceti instabil- 
ity strip. Considering these stars are multi-periodic pulsators and the pulsations 
propagate to the nucleus of the star, they carry information on the structure of the 
star and evolution of the progenitors. The ZZ Cetis discovered till 2003 are mainly 
within 100 pc from the Sun, and probe only the solar vicinity. The recently dis- 
covered ones, and those reported here, may sample a distinct population as they 
were selected mainly perpendicular to the galactic disk and cover a distance up 
to » 400 pc. 
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1. Introduction 
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The ZZ Ceti stars are pulsating white dwar f stars with an atmosphere of pure hydrogen 



McGraw fc Robinson 



1977). They show multi-periodic oscillations with 



periods from 70s to 1500s and fractional amplitudes ranging from 0.4% to 30%. They 
undergo g-mode pulsations caused by the re-7 mechanisms and the development of a sub- 
surface convection zone due to the opacity bump caused by partial ionization of hydrogen 
that starts when the cooling white dwarf reaches effective temperatures around 12 000K. 
The convection zone stores and enhances the heat exchange due to the pulsations. 

The ZZ Ceti class of variable stars is also called DAVs and is the coolest of the three 
known instability strips in the white dwarf cooling sequence: the pulsating PG 1159 



stars, around 200000-65 000 K (Drcizlcr ct al. 1998; 



2004), the DBV, around 25 000-22000 K (Bear 



Naeel fe Werner 2004: 



ichamrj^et^l 



2005a|) and the DAV, around 12270-10850 K l|Bergeron et al 



1999 



Quirion et al 



Castanheira et al 



2004: 



Mukadam et al 



2004). 



The ZZ Ceti class presents gradations between the two extremes: the hot DAVs 
(hDAVs), close to the blue edge of the instability strip, have sinusoidal light-curves with 
low amplitude (<2%) and short periods (< 300s). The cool DAVs (cDAVs), close to 
the red edge of the instability strip, show large amplitude (< 30%) long period pul- 
sations (< 1500 s), non- sinusoidal light curves because they are distort e d by the ex- 



tended convection zone (Brickhil 



1992: 



Wu 



2001 



Ising fe Koester 



2001 



Montgomery 



2004J). Another important factor in shaping the light curve and def ining which periods 



Montgomery fc Winget 



1999; 



Metcalfe et al 



2004; 



ion of the core 




Winget et al. 


1997 


Kanaan et al. 


2005 


), which for the 



high mass (above 1 Mq) white dwarf stars occurs while the star is within the ZZ Ceti 
instability strip, or before it reaches the strip, depending on the mass and the core com- 
position. Pulsations cannot propagate inside a crystallized core, distorting the period 
distribution and decreasing the pulsation amplitudes. 



siars (MrCook fc SioXo0 3 1 . but 



2005; 



There are to date 89 known non-interacting ZZ Ceti stars ( Mullallv etal 

: 11 — 1 

Castanheira ct al. 2005b) among m ore than 5400 spec troscopicall y identified white dwarf 



McGraw 



Coxl ((1980) already indicated they 



1 19771) and I 

are the most common variable star known. Because they are intrinsically faint, My ~ 12, 
the published ones till 2003 are mainly within 100 pc from the Sun, and probe only the 
solar vicinity. The recently discovered ones, and those reported here, may sample a dis- 
tinct population as they were selected mainly perpendicular to the galactic disk and 



Send offprint requests to: kepler@if.ufrgs.br 
* Based on observations at the Southern Astrophysical Research telescope, a collaboration 

between CNPq-Brazil, NO AO, UNC and MSU. 
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co ver a distance up to ~ 400 pc, and the thin disk scale height extends to « 300 pc 



I Maiewski &: Sieeel 



2002) 



White dwarf star s are the end point s of evolution of stars in the main sequence up to 
Wcideman f all stars. Taking into account 

Kepleifl984 ) are global pulsations, with 



around 10.5 Mq (e.g. 



the observed non-radial <?-mode pulsations (e.g. 



each pulsation mode constraining the stellar structure in a different way, we can use 
the pulsations to 



Kepler et al 



2003; 



Hi99oll2nn.4 



untangle the s tructure of the whole star llWinget et al 



Metcalfe 



2003) and even their rates of evolution ( Winget et al 



1991, 



1994: 



1985 



Costa et akin 9991 



Kepl er et al.ll2 000. 



2005 



Mukadam et alJl2003h . These measured 



evolutionary rates have been used to calculate the age of the c oolest known white d warf 



stars, allowing an es timative of the age of the galactic disk IjWinget et al 



of a globular cluster IjHansen et al 



1987) and 



2002). Even more important, pulsating white dwarf 



stars are ex cellent laboratories for testing high energ y and high density physics, suc h 



as neutrino dWinget et al 



crystallization IjWinget et al 



2004 ) and axion emission ( Corsico et al 




2001; 



Kepler 



2004), 



2003), and even an estimation 



Montgo mery et al. 

I.1 I2OO3I . a rate important from early Universe 



of C(a,7)0 reaction rate (Mctcalfi 
composition to supernova explosions, and which determines the size of the C/O core 
of most white dwarf stars. Crystallization, axion emission and cooling rates are mainly 
determined from the study of the DAVs. From evolutionary models, white dwarf stars 
with masses below 0.45 Mm should have He cores, and those above 1.1 Mq should 
have O-Ne-Mg cores (e.g. IWeidemannll2003l) . Another important use of pulsations is to 



use the light travel time variations measurable by the phase changes in the pulsation 
modes to detect planetary companions to the white dwarf st ars. As most planets will 
survive post-main sequence m ass loss to the white dwarf phase ( Duncan fc Lissauerll998 ; 
Mugrauer fc NeuhauseJl2005 ) , we can use the same technique used to study companions 
in pulsars to detect even planets smaller than the Doppler technique can, complementing 
their search space. But planet searches around white dwarf stars require very stable 



2003), and only a small sample 



pulsations, like those found in hot DAVs IjWinget et al 
of them is known to date. 

We are therefore involved in a program to find a significant number of pulsating white 
dwarf stars, to study their structure through asteroseismology, measure their evolutionary 
rates, and look for planets orbiting them. 



2. Candidate Selection 



Survey and fitted to Koester's model atmospheres (Klcinman et al 



acquired by the 


Sloan Digital Sky 


Klcinman et al. 


2004) are good se- 



lection criteria to choose candidates for time series photometric searches of ZZ Ceti stars 



(Mukadam et al 



2004). The SDSS spectra have in general SNR~ 30 for g ~ 18 and we 
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fitted Koester' s spectra models from 3800 A to 7200 A (IKleinman et al.lbo04l) . Unlike the 



fitting done by 



Bergeron et al 



1995, 200 



Koester & Holbcrs (2001), which only fit the 



line profiles and not the continuum, we used the whole spectra from 3800 to 7200A. The 
long wavelength baseline, coupled with the SDSS photometric data, and a low order mul- 
tiplicative polynomial to allow for small flux calibration uncertainties, result in accurate 
T c ff. The selection of this limited wav elength range is to incre a se the weight o f the re gion 
with lines, which are log q dependent. 



Castanheira et al 



Mukadam et al 



Mullallv et al 



and 



I 2005bj) show that we can attain 90% probability of variability if we 



constrain the search range to 11800K > T e g > 10850 K. 



3. Observations 

We used the SOAR Optical Imager, a mosaic of two EEV 2048x4096 CCDs, thinned 
and back illuminated, with an efficiency around 73% at 4000 A, to acquire time series 
photometry. It covers a field of 5.26' x 5.26' on the sky, on the bent cassegrain port of the 
4.1m SOAR telescope. We observed from March to May 2005, when the telescope and 
imager were still under commissioning, even lacking baffle tubes and therefore with an 
increased background. We observed in fast readout mode, with the CCDs binned 2x2, 
which resulted in a pixel scale of 0.154 arcsec/pixel and a readout +write time of 10.2 s. 
The exposure times ranged from 20 s to 40 s, longer than the overhead but still keeping 
the Nyquist frequency in range with the shortest pulsation periods detected to date. 
The data was bias subtracted and flat fielded before we obtained differential photometry 
through weighted apertures around 2 FWHM (full width at half maximum of the seeing 
disk), chosen for highest SNR. We observed each star twice for around 2h each time 
All observations were obtained with a Johnson B filter, considering 



Robinson et al 



1 1982) show the pulsation amplitude increases to the blue, and to minimize the back- 
ground. 

Table [3] list the new variables and their effective temperat ures obtained by fittin g 



the optical spectra to Detlev Koester's model atmospheres, as in 



Kleinman et al 



2004). 



The first pulsator we observed was also observed with Argos l)Nather fc Mukadam 2004) 
at McDonald Observatory 2.1 m telescope, to confirm all the observed periodicities we 
detected, and check validity of the whole observing system (the telescope, the instrument 
and the software). 



4. Results 

In Fig. ^ we show the light curves on the left panels and the Fourier transform of 
them in the right panels. In Table 0] we list all runs obtained for each star and the main 
periodicities detected, i.e., those with a false alarm probability smaller than 1%. {A) is the 
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Table 1. New ZZ Cetis 



SDSS spSpec 
MJD-Plate-Fiber 




Name 


f 


■ 


T cB (K) 


log g 


Main Periodicity a 


52642-1185- 


085 


WD 


0825+0329 


17 


48 


11801+ 105 


8.33+ 0.044 


481s( 


34.5mma 


52650-1188- 


191 


WD 


0843+0431 


17 


93 


11250+ 63 


8.18+0.044 


373s@ 


10.43mma 


52670-1190- 


322 


WD 


0851+0605 


17 


08 


11306 ± 48 


8.11 ± 0.029 


326s<£ 


i22.4mma 


52238-0566- 


031 


WD 


0911+0310 


18 


41 


11634+ 126 


8.11+ 0.084 


347sC 


S17.4mma 


52976-1301- 


445 


WD 


0917+0926 


18 


09 


11341+ 64 


8.15+ 0.044 


289s<£ 


i!16.1mma 


51900-0278- 


367 


WD 


1106+0115 


18 


37 


10990+ 62 


8.09+ 0,049 


822s<£ 


«12.2mma 


52672-1230- 


188 


WD 


1216+0922 


18 


56 


11293+ 109 


8.29 ± 0.078 


823s(£ 


i!45.2mma 


52000-0288- 


412 


WD 


1218+0042 


18 


71 


11123+ 93 


8.16 ± 0.068 


258s 


Q!16mma 


52313-0333- 


077 


WD 


1222-0243 


16 


74 


11398+ 44 


8.35+ 0.026 


396s(£ 


«22.0mma 


52026-0523- 


186 


WD 


1255+0211 


19 


09 


11385+ 154 


8.16 ± 0.106 


897s<£ 


*31.7mma 


51689-0293- 


603 


WD 


1301+0107 


16 


30 


11099+ 34 


8.11 ± 0.023 


879s 


il3mma 


51692-0339- 


629 


WD 


1310-0159 


17 


67 


10992+ 65 


7.92 ± 0.049 


280s(£ 


i!6.56mma 


51955-0298- 


608 


WD 


1337+0104 


18 


57 


11533+ 156 


8.55+ 0.085 


797s(£ 


il0.2mma 


52045-0582- 


551 


WD 


1408+0445 


17 


93 


10938+ 64 


8.06+ 0.044 


849s<£ 


«24.3mma 



a mma is milli-modulation amplitude, corresponding to 1000 x AF/F, where F is the 
measured flux. The MJD-Plate-Fiber are the parameters necessary to access the spectra at 
http:/ /das. sdss.org 



square root of the averag e power, and is an estimate of the noise Ipchwarzenberg-Czernv 



1991 


1999 


Kepler 


1993) 



Mukadam et al 



(2004), 



We also observed the g=18.71 variable discovered by 
WD 1502-0001, SDSS J150207.02-000147.1. Its SDSS spectrum spSpec-51616-0310-206 
fits Toff = 11200 ± 117, and \ogg = 8.00 ± 0.079. with auto23. Its spSpec-51990-0310- 
229 spectrum fits T eS = 11116 ± 096, and logg = 8.18 ± 0.06, with a uto21. Auto23 is 



the ve rsion o f the spectra fitt i ng pro gram and calibration published by 



2004), while 



Mukadam et al 



12004) published values are older auto21. 



Kleinman et al 



Mukadam et al 



2004) measured periodicities at: 687.5s@12.0mma, 629.5s@32.6mma, 581.9s@ll.l mma, 



418.2s@14.9mma, and 313.6s@13.1 mma, classifying it as a cDAV. 

We have also found one star not observed to vary (NOV), WD1359-0034, with a 
detection limit o f 3(A) = 2 mma. Its T c ff = 10640+32 K is outside the main strip found by 



Mukadam et al 



1 2004) , but hotter than their coolest variable. We also confirmed one of 
the NOVs within the instability strip reported by Mukadam et al. 2004, WD 1432+0146, 
with T c ff = 11255 + 73K and log 5 = 8.05 + 0.05, at the detection limit of 3(A) = 4 mma. 
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Table 2. Periodicities detected in the light curves. 



Name 


Mean noise 


Period® Amplitude 


Date of Obs. 


Length 


WD flS^+O^Q 


/ /4 \ — 1 fimm^ 

\/l / _L . VJllilllCL 


481«(Q4 5mma — 9 8/4\ 

ttu X D ' Ni . tjllllllclj Zj . O \ -T7L / 


1 Mar 05 


2h 






303q®3 Sitima 

OVJOO v=itJ . OlllllltX 








(A) — 3 37mma 

\ / O - ' } 1 111 111 CX 


664s@10 7Qmma— 3 2(A) 

UU^tO^ll'- 1 i_/llllllcL O • Zj \ JT. / 


1 1 Mar 05 

_L 1_ IVlcLl \J 'J 


lh 






334s@6 Qmma 

UOtiu vJ ■ 111111CL 








(A)=2.13mma 


644s@12.0mma= 5.63(A) 


13 Aor 05 


2.9h 






704s@6 Omma 










826s@5 3mma 

(J Zj 1J KJ .iJllllllCX 






wd 0843+0431 


I A\ — 3 75mma 


373«(Q)10 4%nma — 9 78//U 

O 1 OoVii 1U .lOlllIllCt Zj . I OV-Tl-/ 


21 Mar 05 

_■ 1_ IvlOil KJ'.J 


2h 




/4\ — 2 09mma 

\/l / Z. . L/Zj111111CL 


10405(0)11 4mma — ^ 64/4\ 

lUtt-Jo ii _L _L . ^±llllllclj tj . VJt: \ Ji 1 


Q Anr O^i 


2h 




//U — 1 Q2mma 

\ / _1_ • i_/ijllllllcL 


1085s@7 42mma— 3 86( , /U 

A-\J\J" J O ii 1 .^tZjllllllCL »J . O VJ \ -/A. / 


1 1 Anr 05 


2.6h 


wn 0851 +0605 


/ A\ — 4 1 9mma 

\/l / 4: . ± ZjllllllcX 


396s<3)99 4mma — 5 4/4\ 

OZj VJo v=iZjZj . ^tlllllltX <J . \ -ii / 


21 Mar 05 


2h 


wd riQi i +031 n 


/ 4\ — 4 ^mma 

\ / j: * OllllllcL 


347q(Q1 7 4mma 

O^t 1 Ov^; _1_ 1 *t:111111cV 


1 Mar 05 


1.9h 






7^7^(0)1 4mma 

( i.J 1 o^ii _1_\J .^tllllllCL 










HSSstail 2 3mma 








/ /4 \ — ^ 1 m m a 

\/l / tj . ± ill ill CL 


3 c i3q(Q)9fi Omma— ^ 1( A\ 


1 1 Mar 05 


1.9h 






17fic<(S)11 Imma 








McD (.4) =4.4mma 


352s@27 7mma— 6 2(A) 


11 Mar 05 


1.6h 






420s(Q1 2 (Smma 






WD 0017-1-0026 


/ A\ — 9 7mma 


98Q«®1fi 1mma— fi/4\ 


1 A A/Tar 0^ 

_L4: IVlcLl U>J 


2.9h 






2501(8)1 2mma 










91 9s(3)S Omma 








( A\ — S Omnia 

\/l / — O . Ullllllcl 


9SS«(S)14 Omma— 4 fi/4\ 


1 ^ Anr 0^ 


2h 


wd i 106+01 is 


(A)— S Omnia 

\ /l / tj . WllllllCX 


899siS)19 2mma— 4 1 (A) 


1 5 Mar 05 


2h 






0S0s(S)10 Imma 








( A\— S Qmma 

\/l / — O . c/llllllct 


037«(S)1 1 Imma— 9S5/4\ 


21 Mar 05 


2h 






71 Qq(Q)7 7m m a 

1 1 . 1 111111CI 










220s@S Omma 

^ ^ \ / O O ■ -Lf llllllct 








all dataZ/H — 9 Fvmrna 


073s(9)1 Smma 

tj 1 UotS _L\J . LJllllllcL 










S49s(3)Q 4mma 

O^tZiOvrs C7 .^tllllllCl 






WD 1 21 6+0022 


( A\— 1 7mma 


S93si0)45 2mma— 4 2M\ 

UtiUov^^lu ■ ^111111<X t: ■ £j \ / 


23 Mar 05 

Zj \J IVlcLl KJ *J 


1.9h 






400s@30 Imma 








M\— 7 Q2mma 

\J1./ i - - / — J.i.L J.llt I 


840s(Bi42 Omma— 5 3M\ 

UTIVO^i^tZl . Wllllllcl tj .O \ il / 


1 3 Anr 05 


2.5h 






570s(S)24 fimma 










626=1(0)21 (Smma 

V^/ Zj V^/o _L <\ w fllllllcL 










967s@20.5mma 






WD 1218+0042 


(A)=4.3mma 


258s@16mma=3. 75(A) 


6 Apr 05 


3.4h 






175s@10.0mma 










100s@ll.0mma 








(A)=2.55mma 


259s@8.2mma=3.2(A) 


15 Apr 05 


2h 






152s@5.1mma 
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Table 3. Periodicities detected in the light curves (cont.). 
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Name 


Mean noise 


Period® Amplitude 


Date of Obs. 


Length 


WD 1222-0243 


(4)=3.1mma 


396s@22.0mma= 7.1(A) 


10 Mar 05 


2h 






198s@7.3mma 








(A)=2.46mma 


198s@6.67mma= 2.7(A) 


11 Mar 05 


Hi 


WD 1255+0211 


(A)=4.84mma 


897s@31.7mma=6. 55(A) 


9 Apr 05 


3.4h 






1002s@21.7mma 










812s@16.4mma 






WD 1301+0107 


(A)=3.8mma 


879s@13mma= 3.4(A) 


11 Mar 05 


1.9h 




(i)=7.8mma 


901s@24mma= 3(A) 


10 Mar 05 


0.6h 




(A)=4.4mma 


870s@22.3mma= 5(A) 


14 Mar 05 


2.9h 




all data (A)=3.0mma 


882s@17.6mma 










628s@15.2mma 






WD 1310-0159 


(A)=2.92mma 


280s@6.56mma=2. 25(A) 


23 Mar 05 


1.25h 






310s@6.4mma 








(A)=2.81mma 


349.6s@17.6mma=6.3(A) 


13 Apr 


2.5h 






280s@9.2mma 






WD 1337+0104 


(A)=3.7mma 


797s@10.2mma= 2.75(A) 


15 Mar 05 


2.5h 




(A)=5.75mma 


735s@18.49mma= 3.2(A) 


i~l -1 TV IT r\ f 

21 Mar 05 


1.6h 




all data(A)=3. Omnia 


715s@10.0mma 






WD 1408+0445 


(A)=4.34mma 


849s@24.3mma=(A) 


15 Apr 


2h 






1038s@12.0mma 










764s@ll.lmma 






WD 1359-0034 


(A)=0.69mma 


NOV2 


6 May 


2.2h 




(A)=2.4mma 




8 May 


l.lh 


WD 1432+0146 1 


(A) = 1.25mma 


NOV3 


11 Apr 


3.6h 



5. Conclusions 



We fit the optical spectra acquired by SDSS with Koester's model atmospheres, deriving 
the effective temperature of the DAs. Selecting to observe with time serie s photometry 



Mukadam ct al 



(2004), we detected 



those inside the ZZ Ceti instability strip derived by 
fourteen new ZZ Cetis, i.e., hydrogen atmosphere pulsating white dwarf stars, in the 
range 11850 K > T cff > 10850 K. 

We do note however that there are 109 stars for which DR3 SDSS have multiple 
spectra, just from 13000 K > T e g > 10000 K, and the fitting results show that the 
mean uncertainties is — 300 K, and eri ogs — 0.21 dex, for the same object but 

different observations. This is larger than the internal uncertainty of the fits, but in 
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New DAVs found with SOAR 



E 1 1 1 1 1 = 
1 — 1 1 1 1 hrH 1 1 1 1 1 1 1 1 1 — ! 


=H Kr 


WDO043+O4O1 ~ 

H — -f — ^tr^yy-y^y^y YN^j^jV^jyvY^ 


j — I — i — i — — i — i — I — — i — i — i — — i — i — i — — i — i — i — — j 




WD 1408+0445 




L /Y/vVvV (V 


WD1310-0159 -= 


= — i — i — i — — i — i — i — — i — i — i — — i — i — i — — i — i — i — — : 
E — I — l — l — — l — l — I — — l — l — l — — l — l — I — — l — l — l — — E 


= Ij , WD1Z55+0211 = 








WD1301+0107 J 


i — I — i — i — — i — i — H--J — i — i — i — — i — i — i — — i — i — i — — = 

p^r *>i -v.- 5»yv «*7 % « 


~ 1 1 1 1 


WD0911+0310 ; 

/ 1 v |/ wyv | — ( — \j — ( — j^i^Y^r ^ 


| — I — h— i — — i — I — I — — I — i — i — [ — i — i — i — — I — I — i — — j 


CI. 


WD1222-0Z43 J 


! — 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 1 






WD0825+0329 : 


^ — I — I — — I — I — I — — I — I — I — — I — I — I — — I — I — I — — 1 
I — 1 1 1 1 1 1 1 1 1 1 h*H h—\ 1 — j 




1 WD0917+0926 : 

" W WUwv/vv^ 


j — I — i — i — — i — i — I — — i — KtH — — i — i — i — — i — i — i — — j 


1 U 1 1 1 1 1 1 1 1 1 1 II 1 1 

— i. 1 1 1 WDiioe+ons - 


• • ™\ •«» •«%••« 

i i i i i = 


-A WD1337+01G4 


: 1 1 1 1 1 1 1 P I* | t 1 1 1 1 1 1 § 

i — 1 1 1 1 1 1 — 1 1 1 — 1 1 1 1 1 1 — — i 


\ /VWl/VVwU 


WD1218+0042 : 


\ 1 1 1 1 1 1 1 1 — H 1 T 1 1 1 1 : 


K-H 1 E- 


. WD0851+0605 J 


I*v\a^^»^*v»^ 




WD1216 + 0932 _| 



10 

5 





15 
10 

5 




15 
10 

5 




10 c 

5 a 



15 
10 

5 




15 
10 
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o.e 
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4000 6000 
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4000 6000 
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Fig. 1. Light curves (left panels) and Fourier transforms (right panels) for the new ZZ 
Cetis. mma is milli-modulation amplitude, corresponding to 1000 x AF/F, where F is 
the measured flux. 



Kleinman et al 



(2004). 



general within 3a of each other - and mostly within 1 or 2a, as in 
The uncertainties cover a substantial fraction of the instability strip. To really study the 
purity of the instability strip we need to reduce the uncertainties to less than 200K, but 
we must take into account the fact that the large amplitude pulsators at the red edge 



have temperature excursions of around 500 K during one pulsation cycle (Robinson et al 



19821) . As the SDSS spectra on average are 3 x 900s exposures per observation, it is unlikely 
that pulsations are causing the 300K differences. 
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Table 4. Periodicities detected for the Variable WD 1502-0001. 
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Name 


Mean noise 


Period® Amplitude 


Date of Obs. 


Length 


WD 1502-0001 


(i)=71mma 


603s@28.1mma=3.95{,4) 


6 Apr 05 


2.7h 






658s@27.0mma 










415s@16.7mma 










141s@14.9mma 








<yl)=7.1mma 


604s@35.5mma=5(A) 


10 Apr 05 


2h 






424s@16.4mma 








Fig. 2. Plot of effective temperature and log g of the observed objects in this paper (filled 
symbols) and SDSS DAVs in general (crosses). The two NOV studied in this paper are 
represented by open rectangles. 
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